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Abstract 

We study gluon production in high energy proton-nucleus collisions in the 
semi-classical hamework of the Color Glass Condensate. We develop a 
general formalism to compute gluon fields in covariant gauge to lowest or- 
der in the classical field of the proton and to all orders in the classical field 
of the nucleus. The use of the covariant gauge makes the diagrammatic 
interpretation of the solution more transparnt. fcx-factorization holds to 
this order for gluon production ~ Our results for the gluon distribution are 
equivalent to the prior diagrammatic analysis of Kovchegov and Mueller. 
We also show that these results are equivalent to the computation of gluon 
production by Dumitru and McLerran in the Fock-Schwinger gauge. We 
demonstrate how the Cronin effect arises in this approach, and examine 
its behavior in the two extreme limits of a) no small-x quantum evolution, 
and b) fully saturated quantum evolution. In both cases, the formalism 
reduces to Glauber's formalism of multiple scatterings. We comment on 
the possible implications of this study for the interpretation of the recent 
results on Deuteron-Gold collisions at the Relativistic Heavy Ion Collider 
(RHIC). 
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1 Introduction 



The physics of high energy proton/deuteron-nucleus collisions has acquired new 
vigor with on-going experiments on Deuteron-Gold collisions at center of mass 
energies per nucleon of \/s = 200 GeV being conducted at Brookhaven's Rela- 
tivistic Heavy Ion Collider (RHIC). First results from these experiments have 
already been submitted for publication and presented at conferences [1-4]. In 
the near future, CERN's Large Hadron Collider (LHC) will study proton-nucleus 
collisions with center of mass energies per nucleon of \/s = 5.5 TeV [5-7]. When 
one considers this important increase in the energy range for these collisions, 
it is interesting to examine how results from lower energy proton-nucleus colli- 
sions are modified both qualitatively and quantitatively as one goes to higher 
energies. 

At very high energies, or small x values, the relevant parton densities in the 
proton and in nuclei grow very rapidly. If these densities are sufficiently large, 
the parton distributions may saturate, [8-10] leading to a qualitatively different 
behavior of the distributions. Saturation will occur sooner in nuclei than in 
protons because the large number of nucleons give rise to an enhanced parton 
density in the transverse plane by a factor ~ A^^^. Proton- nucleus collisions 
therefore, for a wide energy range, provide an attractive physical environment 
wherein the proton probe may be considered as a dilute parton gas with prop- 
erties that are believed to be well understood while the nucleus exists in a 
novel, saturated high density state. This latter state has been called a Color 
Glass Condensate (CGC) and its distinctive features have been extensively ex- 
plored [11-30]. 

In this work, we will formulate a description of high energy proton-nucleus 
collisions within the CGC framework. In this framework, the proton and the 
nucleus are effectively described as static random sources of color charge on the 
light-cone. The leading contribution to particle production is obtained by calcu- 
lating the classical Yang-Mills field created by these sources and then averaging 
over a random distribution. Thus on is led to solve the Yang-Mills equations 
in the presence of two number densities of color charges, pp{x±) and p^{x±) 
respectively for the proton and the nucleus, localized on the light-cone ^. One 
has 

where 

J: = gS''+S{x-)p,Axi-) + gS^'6{x+)p,Ax^) . (1) 

Operators calculated in this classical background field have to be averaged over 
weight functionals Wp and - representing the distribution of the densities 

■^The validity of this approach is not obvious a priori - the motivations for this approach 
have been discussed extensively in the literature and we will not go into them here [31-37]. 
Here, it will be apparent only a posteriori - from its success in reproducing known results and 
its predictive power. 
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of color charges in the proton and nucleus respectively. One has then 

(O) = I [Dpj,][DpMxl pp]W, [x^,pj 0[pp, pj . (2) 

This averaging procedure is essential in order to restore gauge invariance since 
0[pp,Pj^] is computed in a particular gauge. The arguments and Xq denote 
the scale in x separating the large x static sources from the small x dynamical 
fields. It is also through Wp and W^^ that quantum effects, due to evolution with 
X of the light-cone wave functions of the target and projectile, are incorporated. 
The functionals Wp and obey a Wilson rcnormalization group equation 
(often called the JIMWLK equation [15-27]), which governs their evolution with 
Xq and Xq . This equation reduces to the well-known BFKL equation in the low 
density or large transverse momentum limit. In this paper, we will explicitly 
show that the physical quantities relevant for proton-nucleus collisions can be 
related to correlators of Wilson lines, the quantum evolution of which can be 
studied thanks to the JIMWLK equation. 

The first part of this work, on gluon production in pA collisions, will be 
presented in this paper. The second part, on quark production, will start from 
the solution of the Yang-Mills equations obtained in the present paper and will 
be discussed in an accompanying paper [38]. 

This paper is organized as follows. In section 2, we show how the Yang-Mills 
equations can be tnmcatcd systematically in order to obtain a coupled set of 
equations for the gauge fields to first order in the proton source density pp and to 
all orders in the nuclear source density . This approximation is valid as long 
as the parton densities in the proton remains small: this is of course the case 
at moderate collision energies, but also at very high collision energies albeit in 
a more restricted kinematical domain (not too far away from the fragmentation 
region of the proton). Gauge invariance is preserved, by construction, in this 
truncation. We explore the structure of these truncated equations and introduce 
a convenient diagrammatic expression for the gauge fields. The equations are 
solved explicitly in section 3. We obtain a compact expression for the gauge 
fields in terms of two distinct Wilson lines whose coefficients are shown to be 
simply related to the well known effective Lipatov vertex [39-42]. Classical 
gluon production is discussed in section 4. It is shown explicitly how the results 
can be written in a A;_L-factorized form as the product of two A;_l -dependent 
distributions times a factor proportional to the square of the Lipatov vertex. 
The results obtained in this fashion are compared to the results of Dumitru 
and McLerran [43] in Fock-Schwinger gauge and are shown to be identical. In 
section 5, we show how the Cronin effect arises in simple models with and 
without quantum evolution of the color charge densities. Readers interested in 
this topic alone may skip directly to section 5. Interpreting our results in light 
of the recent data from the BRAHMS experiment at RHIC for both central and 
forward rapidities [1] , we find interesting ramifications for our understanding of 
the physics from both d-Au and Au-Au experiments at RHIC as well as for the 
future Heavy Ion experiments at CERN's Large Hadron Collider (LHC). Various 
technical details are included in appendices A and B. Finally, we comment on the 
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relation of results in our approach to the Glauber multiple scattering formalism. 
Since this is somewhat outside the main thrust of this paper, we discuss it in 
a self-contained fashion in appendix C. We conclude in section 6 and set the 
stage for the results on quark production to be presented in [38] . 

Before we proceed further, we should point out that several results on gluon 
production in this paper are not new. They were first derived, in covariant 
gauge, by Kovchcgov and Mueller [35] and subsequently refined in later pa- 
pers [44,45] . These results were also derived, shortly after the work of Kovchegov 
and Mueller, by Dumitru and McLerran in the Fock-Schwinger or radiation 
gauge [43] . (We should also mention the work of Braun in a somewhat different 
approach [46].) Our derivation, albeit also in covariant gauge, follows a differ- 
ent tack from the Feynman diagram technique of Kovchegov and Mueller. We 
solve the Yang-Mills equations explicitly for the light-cone sources in Eq. 1 . The 
expressions for the gauge fields are derived explicitly and one observes how the 
Lipatov vertex arises in this approach. We also show, for the first time, the 
exact equivalence between the approach of Kovchegov and Mueller with that of 
Dumitru and McLerran. Our results for the gauge fields are new and will be 
useful in deriving the result for quark production, which will be discussed in [38] 
(henceforth referred to as paper II). Our discussion of the Cronin effect is also, in 
part, known from previous work by us and other authors [47-52,45,53-55]. Our 
results in the "super-saturated" limit are new and help clarify the interpretation 
of the RHIC data. The fact that one recovers the Glauber multiple scattering 
picture in the case where the correlations of the color sources are Gaussian is 
well known [56-58] - however the formalism discussed here provides a system- 
atic way to go beyond the limitations of the Glauber model by incorporating 
the non-Gaussian correlations that may arise via quantum evolution. 

2 Yang-Mills equations 

The Color Glass Condensate is a classical effective field theory which describes 
the physics of high energy, semi-hard processes in QCD. In this approach, one 
first solves the classical Yang-Alills equations in the presence of the light-cone 
sources pp and , next one computes the observable of interest in this classical 
field, and finally one averages over the sources pp and p^. For proton- nucleus 
collisions, one requires the solution of the gauge field to lowest order^ in the 
source density Pp of the proton and to all orders in the nuclear source density 
p^ . We shall write down in this section the truncated set of Yang-Mills equations 
to this order. We next show that these equations can be written in a compact 
form and explore their structure and their diagrammatic content. Solutions to 
the Yang-Mills equations will be discussed in the next section. 

Going beyond this approximation and solving the Yang-Mills equations to all orders in 
both sources has only been possible numerically so far [36,37,59-61]. 
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2.1 Generalities and notations 

The classical Yang-Mills equations read*: 

[D^,Fi^^] = r. (3) 

J" is a color current which at lowest order in the sources pp and reads: 

J: = g5''+5{x-)pp,^{x,_) + g5^-6{x+)p,,^{x,_) . (4) 

Pp is the color source describing the proton, moving in the +z direction at the 
speed of light, is the color source describing the nucleus, moving in the oppo- 
site light-cone direction. These sources represent the number densities of color 
charges in the proton and nucleus respectively, which are Lorentz contracted to 
delta functions [31-33] - 5{x^). In general, the current receives higher order 
corrections in Pp,p^ because it must be covariantly conserved: 

[D,,r]=0, (5) 

and the produced gauge field has a feedback on the current itself. Moreover, 
this system of equations is still under-determined because of gauge invariance. 
Here, we choose the Lorenz gauge: 

a^A" = . (6) 

Using this gauge condition, one can bring the Yang-Mills equations to the fol- 
lowing form 

□A" = r + tg[A^„ F^" + 5^ A"] . (7) 

Since the commutator starts at least one order in pp_^ after the gauge field itself, 
this form is appropriate for an expansion in powers of the color sources. 

In this section where we deal with truncations of the Yang-Mills equations, 
it is convenient to introduce the following notations. For any quantity X, let 
us denote by X(m,n) the term of order p™/0" in its expansion in powers of the 
sources. We therefore have: 

+ 00 +00 

^ = XI X] ^im,n) ■ (8) 
m— n— 

Let us also denote by X(o_oo) and -''^(i.oo) the quantity X at order and 1 in pp 
and to all orders in p^, i.e.: 

-|-oo 

-'^(O.oo) = ^ ^(0,n) ) 

n=0 

+ 00 
71=0 

Thus, for pA collisions in a kinematical regime where the proton can be seen as 
a dilute object, we want to determine the component of the gauge field. 

*By convention: D^^ = — igAl^. 
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2.2 Order p° 

Before we proceed to consider the field generated in proton-nucleus collisions, let 
us first discuss the well-known case of the classical field of a nucleus (before the 
collision). If we first rewrite the Yang-Mills and current conservation equations 
to order p^, we have: 

n^rO,cx>) = -/(O.oo) +i5[Ao,oo)^,i^fo';x>) +^''^ro,oo)] , 

^M"^(0,oo) = «fl[^(0,oo)^, ^("0,00)] • (10) 

This system of equations can be solved iteratively to the appropriate order in 
pj^. Indeed, at any given step, the commutators of the r.h.s. always involve the 
gauge field obtained at the previous step. At order , we have simply: 

J>;;^,^=g5>^-5{x+)pM, 

D^(o,i)=-^ro,i)- (11) 

Note that the current conservation is automatically satisfied at this order, 
because has only a — component that does not depend on x~ (hence^ 
9^ J^Q 1) = ^^-^(o 1) =0). The solution of the Yang-Mills equation is also trivial^: 

Al^,^=-g5'^-5{x+)^p,{x^). (12) 

Note that the field strength corresponding to this solution has only two non- 
vanishing components: 

-^(0,1) ^ ""^(o.'i) ^ '^'^(0,1) ■ (1"^) 
At the order , the equations are: 

n^(o,2) = Jio,2) + i5[^(o,i)^, i^('^o:i) + ^''^(0,1)] • (14) 

Given the solution obtained at the first order, the commutator [A(o,i)m' "^{o i)] 

that appears in the current conservation equation vanishes, and the solution is 

"^(0 2) ^ if we impose that the current vanishes in the remote past. The same 

is true of the commutator [A(o,i)^, -F(o\) + d^A'^^ ^^j, so that we have simply 

□A^Q 2) ~ 0. Again, if we impose that the field vanishes in the remote past, we 

= d/dx-, d- = d/dx+. 
^This notation is a convenient abuse of language. One should of course read; 



{x±) = J d'^y±{'^±\:^\y±)pAiy±) 
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have Aj^Q 2) ~ 0- Therefore, the gauge field has no term of order PpP^. These 
arguments can be trivially extended to any order in , so that the solution at 
order p^pp\ is in fact valid to all orders in : 

^fo.oc) = ^fo,!) = -95'-K^^)^^pA^1-) ■ (15) 

In section 3 and following, we will denote the field of the nucleus alone simply 
A^^ instead of ^{^'g . Similarly, the field of the proton alone will be denoted 
instead of A'^^ . 

2.3 Order pp 

We now come to the case of physical interest. At the order p^, the equations 
we need to solve are the following: 

'-'^(1,00) — "^(1,00) 

+zff[Ai.oo)M,^ro:i) +^''^^0,1)] +^5[^(o,i)M,^^ri:oo) +^'^^(1,00)] , 

9t^J{i,oo) = ^5[Ai,oo)M. -^fo,!)] + ^9[-4(o,i)M, J{i.^oo)] ■ (16) 

Let us begin with the equation for the current. We can make the reasonable 
assumption^ that J^^ = 0, i.e. the current J^^^ is purely longitudinal. This 
assumes that the sources that produce this current do not undergo any recoil, 
i.e. that they are eikonal light-cone sources. In addition, we require that the 
component , which is a correction from interactions with the nucleus to 

the current produced by the proton, depends locally^ on pp{x±) (as was 
Since in the r.h.s. of the current conservation equation, the only term which is 
local in pp is the term ig[A(o,i)^i, J({ 00)]' ^® '^^^ ^^^^ ^P^^* equation into 
separate equations for the + and — components of the current, as follows: 

^ '^(1,00) = ^5'[^(o,i)' "^(1,00)] ' 

9^J[i,oo)=i9[A+^oorJ[o,i)]- (17) 

We can then write more explicitly the evolution equation for each component 
of the gauge field A'^^ , using all the information we know from the field at 

the previous order. For we get: 

□4,00) = 4,00) +^5[^ro,i)'^^4,oc)] • (18) 

For the transverse components A^^^ , we have the following equation: 

n^koo) = «fl[S'^ro.l)'^("l,oo)] - ^5[^r0.1)'^'^("l.oo)] + 2i<7[^^o,l)'^^^(l,oo)] • 

(19) 

■^The fact that we can find an exact solution of the Yang-Mills and current conservation 

equations under this assumption will justify it a posteriori. 

^This means that pp{x±) cannot enter in ^i(asx) via (1/Vj^)pp(a!x)- 
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And finally, the equation for A^^ is: 

+^3[^ro,i)'^(too) + ^-'^a,oo)]- (20) 

This last equation can be made more explicit if we replace the strength tensor 
by its expression (to order pp and all orders in p^): 

^(too) = ^^^a.oo) - 5^4,00) - ^5[4,^), A^o,l)] ' (21) 

which leads to: 

+(^5)'[4o,i)'[4i,cx>)'4o,i)]] 

+2i5[^h,oo)' ^'4,1)] + 2i5[^^o,i)' ^+^a,oo)] • (22) 
2.4 Interpretation of the equations 

In order to make the connection between the solution of the Yang-Mills equa- 
tion and the perturbative diagrammatic expansion, it is useful to interpret the 
equations derived in the previous section in terms of Fcynman diagrams. These 
equations can be simplified by first rewriting all the commutators in terms of 
matrices of the adjoint representation, as follows: 

[A, Bl = ir'-'A,B„ = - {T")^,^ A,B„ = [{A ■ T)B]^ . (23) 

We can then exploit this compact notation to rewrite the current conservation 
and Yang-Mills equations of the previous subsection as follows: 

(^"-^ff4o.l)-^)^(toc)=0' 

^+-^a,oo)= ^5(4,00) •^)-^(o,i)' (24) 

and 

(n-^5^ro,i)-^^^)4,oo) = 4,oc)' 

+i5(5M-,^.r)4_^), 
(n-2i5^^o,i)-^^+)4,oo) = "^a,oo) 

-^9(^ro,i) • ^)^"4,oo) + ^5(S-A^o,i) • 2^)4,00) 
-{i9)\A-^o,i) ■ Tf^too) - 2^5(aM^o,i) • T)^(i,oo) (25) 

When writing these equations, we have collected in the l.h.s. all the homoge- 
neous terms, and kept the inhomogeneous terms in the r.h.s. One then sees 
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that the structure of these equations requires that we solve them in a specific 
order: we need to first determine Jji[ then J^^ and finally 

^(1 oo) • ^^^^ because the quantities previously determined enter as source 
terms in the evolution equation for the latter quantities. 

It is in fact easy to understand how the various components of the gauge 
field are coupled in these Yang-Mills equations by studying the structure of the 
gluon propagator in the color field created by the nucleus. Forgetting 

the source terms which are irrelevant in this discussion, we can formally 

write the field A'^^ oo)(^) ^ "time" .t+ in terms of the field A'^^ cx3)(2/) ^■t 
earlier time Xq via the retarded gluon propagator: 

^(i,oo)(^)= / dy-d'y^G,^4x,y)2d+A1,^^^{y). (26) 

y+=—oo 

where the integration is carried out on an hypersurface = const (taken here 
to be infinitely remote in the past). G^" [x, y) is the retarded gluon propagator 
(in the presence of the background field ^■^). Eq. (26) is not obvious and 
is derived in appendix A. In the following discussion, it will be convenient to 
consider the Fourier transform of the propagator, defined as follows: 

= j d''xd*ye">-^e-'P-yG';^-^{x,y) . (27) 

The full propagator can be obtained by multiple insertions of the external 
field A^Q created by the nucleus. There are two possible insertions, correspond- 
ing respectively to a three gluon vertex or to a four gluon vertex, as illustrated 
in figure 1. The free propagator, proportional to g^"^, does not mix the various 



p,a o,b p,a a,b 

A - C A - C A - d 

^(0,1) ^(0,1) ^(0,1) 



Figure 1: TJie three gluon vertex r^^^(g,p, /c) and the four gluon vertex 
^bacd^ {q,P, fci, fe) that can be inserted on the gluon propagator. 

components ^y and A*^^ of the gauge field. This property how- 

ever does not hold for the full propagator. In order to see this, let us consider 
the correction corresponding to the insertion of one external field A^ , via the 
3-gluon vertex: 

5G^^(9,P) = Gl^^iq) (r^r+(9,p, k)A-:i)('')) GO /(p) . (28) 
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The vertex T"^^^ (q,p,k) is given by the usual QCD Fcynman rules in the co- 
variant gauge, and reads (see figure 1 for the direction of the momenta): 

KSHq,P,k)^gr"^[g^^p + q)++g'^+{-q-kr + g''+{k-pr] . (29) 

More explicitly, its components are given by; 

ri,t+(<z,P,fc) = o, 

T^-^iq,P,k) = -g{p + q)-r''% 

Ka'c^i^^P^ k) = r,-++(<?,P, k) = gp+r"^ , 
Tlt+{q,p,k)=T+:+{q,p,k) = 0, 
T^:t{ci,P,k) = -g{q + kyr'^ , 
Ti-J{q,p,k)=gik-pyr'\ 

^ZU<l,P,k) = 2gp+g'^r'^ . (30) 

In these expressions, we have used the fact that k'^ = (because the background 
field j4j^ does not depend on a;^). This means that all the components of SG'^'^ 

are non-zero, except for (5G^* and SG^^ . Looking back at eq. (26), this 

implies that all the transitions A'^^ — > A'^^ are allowed, except for the 
following transitions: 

^(l,oo) ^(l,oo) ' 
"^(l,oo) ^ "^(l,oo) ' 

4,00) -4,00). (31) 

which cannot happen. The insertion of the 4-gluon vertex does not change this 
conclusion. Indeed, the corresponding correction to the propagator is given by: 

SG'/^{q,p) = Go^/{q) (K^:/ A-^^^,A-^^^,) Go^^(p) , (32) 

with a vertex that reads: 

r''P++ = ig'^SP- S"- {T'^T'^ + T'^T'')ab ■ (33) 

Therefore, only the component 5G'^ is non-zero. This corresponds to a tran- 
sition At ^ — > AT^ 

(l,cx>) (l,oo) 

Our discussion thus far goes a long way towards explaining how the equations 
for A^^^ and A^ are nested and what the couplings between the 

various field components are. We note the following salient points: 

• cannot be produced from A^^^ or from A~^^ Therefore, the 

equation for A^ does not involve any other component of the field, 
and its only source term is the term in J^i ^oy Moreover, the term 
in ig{A'^Q ■ T)9+A^ is the term that corresponds to the transition 
^(1 oo) ~^ ^(1 oo) insertion of one power of the background field ^y 
i.e. to the vertex P"' '". 



10 



• can be produced from Since there is no source term J^^ ^y 
the inhomogeneous terms in the equation for A^-^ corresponds to the 
vertex that is responsible for the transition A^^^ ^y i.e. F* '". 
The homogeneous term in this equation corresponds to the transition 

^1 i-e. to the vertex r"+. 

• In addition to the source term associated with the current XT A7. ^ 
can be produced from (vertices F and F '"+) or from 
(vertex r^*+). The homogeneous term in the equation for A'^^ corre- 
sponds to the transition A^^ — > A'^^ ^y i.e. to the vertex F 

One should note however that the equation for A'^^ ^y as written in 
eq. (25), does not quite fit with this interpretation. (Among other things, 
the homogeneous term is twice too large for this interpretation to hold.) 
However, this can be rescued by subtracting from the r.h.s. of this equa- 
tion the quantity ig[A'^Q ^ydfj^A'^^ ^^], which is zero according to the gauge 
condition we have chosen. After this subtraction has been performed, the 
equation for A^^ becomes: 

(n- ^5^^0,1) •^^^)^ri,oo) = -^a,oo) 

-2z5(A^o,i) • ^)^"4,cx>) - ^fl(^"^ro,i)-^)4,oo) +5'(^ro,i) • ^)'4,cx>) 
-i5(A^o,i) • r)S'^'(i,oo) - 2*5(9^^0,1) • ^)^'(i,oo) • (34) 

One now notes that all the terms in the equation correspond to the above 
interpretation. 



2.5 Diagrammatic content of the equations 

At this stage, it is possible to sketch the Feynman diagrams that correspond to 
the terms included into the solution of the Yang-Mills equations at this order. 
The simplest case is the diagram corresponding to the solution at order . Since 
the field is linear in the source p^, we can represent it by the diagram of 

the figure 2. 

The current is obtained by multiple scatterings of the source pp in 

the field A'^^ ^y as illustrated by the diagram on the left of figure 3. 

A typical contribution to A'^^ is obtained by producing a ghion from the 
current ^y and then by letting this gluon propagate through the background 
field A'^Q . This is illustrated by the diagram in the middle of figure 3. The 
field has a source term proportional to multiplied by a vertex 

F' that enables the transition Af, n Al. . There cannot be more than 
one such vertex. Once formed, the field propagates in the background 

field AVq as illustrated by the diagram on the right of figure 3. 
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Figure 2: Diagrammatic representation of the field A^^^ A solid boldface line 
represents one power of , and a dashed boldface line represents one power of 
Pp- 




(l,oo) 



Figure 3: Diagrammatic representation of the current and of the fields 

^(1 oo) ^'^'^ ^(1 oo)- -^^^ circled vertex is the vertex where the transition 



There are several sources for the field . One possible source is a direct 

transition A^^ as illustrated in figures 4 (there are a 3-gluon and a 

4-gluon vertex for this transition). Another contribution to A'^^ is due to the 
transition A^^^ as illustrated in figure 5. Finally, one last 

way to produce is via the source J^^ , which requires the interaction 

of an with a source . Note that this interaction is local in . The 

diagrams contributing to J^^ ^y as well as the corresponding contribution to 
^(1 cx>)' illustrated in figure 6. 



"5"' 

11 lllll ' 11 llAli ■ 



Figure 4: Contributions to A,. . by a direct transition Af. . A,, Left: 
via the 3-point vertex; the boxed vertex is the vertex T '". Right: via the 
4-point vertex; the dotted boxed vertex is the vertex F '"+. 
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11 llllllll 



Figure 5: Contributions to A^^ by a transition 
The triangular vertex is the vertex r~*+ . 



■3~3 



11 ii\j,r., 11 ii 



Figure 6: Contribution to J^-^ and the corresponding contribution to A^^ 

3 Solution of the gauge field equations 

Now that we have written down the truncated set of coupled equations to lowest 
order pp in the source density of the proton and to all orders in the source 
density of the nucleus, and understood the structure of these equations both 
analytically and diagrammatically, we are ready to solve the equations. This we 
will do in the order prescribed in the previous section. Since at this point, there 
is no risk of confusion, we can simplify our notations by removing the indices 
that indicate the orders in pp and p^. From now on, a quantity relative to 
the proton alone simply carries a subscript p, a quantity relative to the nucleus 
alone has a subscript A, and a quantity at order (1, oo) has no subscript at all. 
In particular, we have the following correspondence between the notations of 
the previous section, and those of the rest of the paper: 

^(0,1) ' ^(1,0) ' 

^fl.oo)--^'' . ^(1,00) • (35) 

3.1 Expression for the current J"*" 

In order to determine the source J+, we have to solve the equation: 

{d- - igA- . T)J+ = , (36) 

with the initial condition J+(a;+ = — oo) = J+ = g6{x~)pp{x±). The solution 
can be expressed as: 

J+(a;+,a;~,x_L) = g U{x~^ , -oo; x±) 6{x~)pp{x±) , (37) 
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where 



U{x^ ,xf;x±) = P+exp 



igjj dz+A-{z+,xi_)-T 



(38) 



is an adjoint Wilson line in the background field A'^ (V+ denotes the ordering 
along a;+). Here A~ can be simply expressed in terms of the nuclear source 
density through the expression given in eq. (15). Note that eq. (37) depends 
locally on pp{x±), but non-locally in p^(a;_L) (because it depends on via the 
Weizacker- Williams field A~). 

3.2 Expression for ^4+ 

The next thing to do is to determine the field A'^. It obeys the following 
equation: 

(□ - igA- ■ Td+)A+ = J+ . (39) 

We will go into the derivation of its solution in some detail because several 
subtle issues arise (regarding smearing of sources and the appearance of two 
different Wilson lines in the solution) that will be relevant for the discussion in 
this paper and in [38]. Using eq. (133) in appendix A, we can write the solution 
in terms of the source and of the initial condition at y+ = — oo as follows: 

A+{x) = ld^yG^ix,y)J+iy) 

+ J dy-d'y^G,{x,y)2d+A+{y). (40) 

y+ — — oo 

where {x, y) is the retarded Green's function for the operator □ — igA~ ■ Td^ , 
defined by: 



(□, - igA-{x) ■ Td+)GJx,y) = 5^^\x - y) , 

x+—y+—>0+ Z 

G^{x,y)cxe{x+ -y+) . (41) 



The limit a;+ — > t/+ can be obtained from the free retarded propagator in the 
Feynman gauge^: 

/j4i, -ik-{x-y) fiv 
(2^FT^ = '^Oix--y-)Oix--y-)Siix-yr) . 

(42) 

^When the time difference goes to zero, the structure of the full propagator in 1+ becomes 
similar to that of the free propagator because there is not enough time for interactions with 
the background field to taJje place. 
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When .T+ — oo, the gauge field must vanish because this is before the projec- 
tiles start contributing to the field. Therefore, we can take the second term in 
eq. (40) to be zero, and write simply: 

A+{x)= jd%G^{x,y)J+{y). (43) 

The background field A~ is proportional to S{x^), but some expressions are 
ambiguous unless we regularize this delta function by giving it some small width. 
Therefore, we are going to replace in intermediate steps the 6{x~^) by a regular 
function 5e(a;+) that is non-zero for a;+ e]0,e[, and such that: 



/' 

Jo 



dx+5^{x+) = l. (44) 



Note that the necessity of such a regularization is a peculiarity of the Lorenz 
gauge. Other oddities of this gauge will appear in the course of our calculation, 
which will naturally disappear from all physical quantities. Along with this 
regularization, we break down eq. (43) into three contributions, according to 
the value of y+: 

A+{x) = f dy+ ( dy-dy^G^{x,y)J+{y) 

J — OO J 

+ ^ dy+ J dy-dy^G^{x,y)J+{y) 

+ J^^ dy+ J dy-dy^G^{x,y)J+{y) . (45) 

We have assumed in this decomposition that a;+ > e. If this were not the case, 
some of the terms may be dropped. 

Let us first discuss the expressions for the current in eq. (45). The expressions 
for the propagator will be discussed subsequently. In the first term, (y+ < 0), 
we can replace the current J+ by J+. In the intermediate term (0 < t/+ < e), 
we have^°: 

J+iy) = g U{y+,0; yj 5{y-)p^{y^) . (46) 

It is clear that the U which appears in this formula corresponds to a color 
rotation of the incoming sources due to interactions with the nucleus. In the 
third term (y"*" > e), we have: 

J+{y)=gU{e,0;y^)6{y-)pp{y^) . (47) 

We now consider the expressions for the propagators in the three terms in 
eq. (45). First we note the following convolution property that can be obtained 
by applying twice eq. (26): 

GJx,y)= J dz-d^z^G,{x,z)2d+G,{z,y), (48) 

z+=const 



^"^Since the support of A is < < e, we have U{x~^, — oo; x±) = U{x~^,0; x±). 



15 



valid for j/+ < z"*" < (see eq. (48)). This may be used in the first term of 
eq. (45), where x~^ > e and < 0, in order to spUt the propagator into three 
pieces: 

G^{x,y)= J dv-d^v^ J dw-d^w^Gl{x,v)2d+G^{v,w)2d+Gl{w,y) , 

(49) 

where we have used the fact that the propagator is a free propagator when its 
two endpoints are both in [— oo, 0] or in [e, +oo]. Similarly, for the second term 
of eq. (45) (a;+ > e and < t/+ < e), we split the propagator into two pieces: 



G'fl(a;,2/) 



J dv-d^v^Gl{x,v)2d+G^{v,y) , 



(50) 



while in the third term (a;+, y+ > e) we simply have: 

G,{x,y)=Gl{x,y). 



(51) 



We see that the only non-trivial propagator we need is the propagator connecting 
two points whose "times" x+ and lie in the range [0,e]. In this case, the 
retarded propagator simply reads: 



GR{x,y) = -e{x -y )d{x± - y_^)V{x+ ,y+;y_^) 



where: 



V{x+,y+;yj_) = P+ exp 



J dz+A-{z+,y^)-T 



(52) 



(53) 



This path ordered exponential is distinguished from the path ordered exponen- 
tial U that we encountered previously by the factor of 1/2 that appears in the 
exponent. Notice further that 



'2d;^G^{x,y) = 6{x -y )6{xjl - yj_)V{x+ ,y+;yj_) . 



(54) 



With the expressions for the currents and the propagators in the three dis- 
tinct kincmatical regions, we can now put together our expression for the gauge 
field in eq. (45). We first use the expression in eq. (54) for the propagator and 
the simple form for the current in the region y+ < to write the contribution 
to from this region as 



9 j dy+ j d^y^ j dv-d^v^Gl{x,v)V{v^)2d+Gl{v,y)pp{y^) , (55) 



y-=0 



v+=0 



where we use a simpler notation for Wilson lines that run over the whole longi- 
tudinal extent of the nucleus: 



V{x±) = V{e,0;x±) = V(-hoo, -oo; a;_L) 



(56) 
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Wherever possible, we have taken the limit e ^ 0. Note that for the 1, i.e. the 
identity term, in the V path ordered exponential, one can trivially convolute 
the two free propagators and obtain: 



9 J dy+ J d^y^Gl{x,y)pp{y^) , (57) 

-oo y~=0 

which is nothing else but a piece of the field of the proton. 

The contribution from the domain j/+ e [0, e] goes to zero when e — > 0. The 
only other contribution to is the one coming from the range > e: 

+ 00 

9 j dy+ j d'y^Gl{x,y)U{y^)pp{y^) . (58) 

y-=0 

The 1 in the ordered exponential U gives the second half of A+. 

We can now combine the above results to write the complete solution for 
A'^{x) in the region where a;+ > 0: 



A+{x) = gJdy+Jd^y^Jdv-d^v^Gl{x,v)V{v^)2d+Gl{v,y)pj,{y^) 

— OO y-=0 v+=0 

+ 00 

+ g jdy+ j d^y^Gl{x,y)U{yMVA.) ■ (59) 

y-=0 

Physically, this result made of only two terms suggests the following. In a col- 
lision at very high energy, the proton must emit the gluon either before or after 
hitting the nucleus. Emitting the gluon during the collision with the nucleus 
is very unlikely at high energy because of Lorentz contraction of longitudinal 
distances. (Such a process is suppressed by at least one power of the inverse 
collision energy.) Note however that, rather surprisingly, both path ordered ex- 
ponentials U and V contribute to the expression of the gauge field. Naively, we 
would have believed only the contribution from the color rotation of the sources, 
the J7's would have been relevant; the V contribution, from the particular form 
of the propagator in this gauge, is unexpected. We will see that our intuitive 
expectations will be restored when we compute the gluon distribution 

In this derivation of . we rightly assumed throughout that > e. How- 
ever, for later use in the calculation of and A~ , we need also the value of 
for x+ e [0, e]. Only the term where the gluon is emitted before the 
interaction of the proton with the nucleus can contribute here, and the last free 
propagator does not appear. We obtain, 



A+{0<x+ <e)=g J dy+ J d^y^Vix+,0;x^)Gl{x,y)pp{y^). (60) 



Albeit, further surprises await when we consider quark production [38]. 
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It is also useful to compute the Fourier transform^^ of the field A'^{x), for 
which we obtain^"^ 



A+{q)=A+{q) 



+ 



[q- + ie){q^ + iq+e) 



Jg2q+_ 
q^ + iq'^e 



(fk 



U{y±) - 1 pp{y±) 



(61) 



This expression can also be put in the more compact form 
A+{q) = A^iq) 



+ 
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q^ + ig+e 



[U{k21_) - {2'Kf5{k21_)\ 



+2g+[F(fe2±)-(27r)25(fe2±)] 



Pi{kii_) 
klj_ 



(62) 



with the following notations: 



k2± = q± 


~ ki± , 










y(fe^) ^ ^ 




■^-V{x^) , 


Ppik±) = 







(63) 

When we expand the path ordered exponentials to first order in A~, we 
recover the solution by Kovchegov and Rischke [33] for the gauge field to lowest 
order in both sources: 



9'^(i,i)(9) = -.9' 



(Pki 



P^(fc2±) -T Pp{ki_L) 



(64) 



Note that it only thanks to the fact that the Wilson line V has a factor 1/2 in 
its exponent that this lowest order result is correctly reproduced. 

3.3 Expression for 

Since A'(a;+ = —00) = 0, the component A* of the field can be expressed as: 

A'{x) = j d\D^{x,v)J'{v) , (65) 

^^We do not use a special symbol for Fourier transforms, since the context makes obvious 
what object we are dealing with. 

l^When computing this Fourier transform, one must remember that eq. (59) is not valid for 
a;+ < 0. In the region a;+ < 0, we simply have A+(a;) = Ap(x). 
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where is the retarded Green's function of the operator □ — 2igA^ ■ T9+: 
- 2igA- {x) ■ Td+)D^ {x,y) = ,5 W {x - y) , (66) 
and where we denote 

J\x) = ig[{d'A-{x) ■ T)A+{x) - {A-{x) ■ T)d'A+{x)] . (67) 

Note that J*(x) is non-zero only for x+ G [0, e] because of the proportionaHty 
to the field A^{x). The propagator is almost the same as the only 
difference is that the ordered exponential V is replaced by a C/ (with the same 
arguments) in the region of interaction with A~. It is easy to check that the 
propagator D^ {x,v) connecting a point in the interaction region (?;+ e [0, e]) 
and a point after the interaction region (x+ > e) is given by: 

D^{x,v) = Gl{x,v)U{e,v+;v^) . (68) 

We need also eq. (60), which we can rewrite as follows: 



A+{x) = gV{x+,0;x^) J dy+ J (fy^Gl{x,y)pp{y^) 

— oo y~=0 

= V{x+,0;xA_)A+{x) . (69) 
We can now put together all the ingredients in order to obtain, 



x[{dlA-{v)-T)V{v+,0:v^)A;{v) 
-{A-{v)-T)dl{V{v+,0;v^)A;{v))] . (70) 
The Fourier transform of this expression is given by 

e 

AUq) = —^±—fdv+dv-(fv^e'''-''U(€,v+;v_L) 
r + J 



x[{dlA-Jv)-T)V{v+,0;v^)A+{v) 
-{A^iv)-T)dl(y{v+.0:v^A^{v))] . (71) 

At this stage, in order to further simplify our result, we must find a way to deal 
with the convolution of a,U and a V. Using the results justified in appendix B, 
we can write this expression as 

x{{V{v^)-U{v^))dlA;{v) + {dlV{v^))A;{v)] . (72) 



i^The factor of 2 in eq. (66) is what distinguishes Dj^ from the propagator in the 
previous subsection - see eq. (41). 
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A final simplification occurs after writing the terms on the r.h.s in terms of their 
Fourier components, 



A\q) = 



2ig 



+ iq+e 



(Pki 



'^{k{{V{k,^)-U{k,^)) + kW{k,^)Y-^^ . (73) 



Expanding the time-ordered exponentials to first order in the background field 
A~ , we recover the lowest order solution of Kovchegov and Rischke [33] for the 
transverse components of the gauge field at that order, 



«'^ii,i)(9) = -5' 



^2 ^1 



PA{k2±) ■ T pp{ki±) 



A-2 



fr2 



(74) 



3.4 Expression for A 

One way of determining A" would be to solve its equation of motion with the 
appropriate source terms. Since there are several source terms, this procedure 
is somewhat lengthy. It is much easier instead to obtain the final component 
of the field A^^ from the gauge condition. In Fourier space, the latter reads as 
follows: 

q-A+iq) + q+A-{q)-q^-A^{q)=0. (75) 
This further implies that 



A-{q) = — [q^-A^iq)-q-A+iq)] . 



(76) 



Substituting our results for and A^ from the previous two subsections, 
we obtain. 



A-{q) = 



q^ + iq+ej {2tt) 



d^kn 



[U{k2±) - {2nf6{k2±) 



+ (2^ - 2g-) [V{k2^) - {2.rS{k2^)] I . (77) 

At order in p^, this expression is zero as expected. When we expand the 
time-ordered exponentials to first order in this formula, we again recover the 
known lowest order result [33], 



(PkiA 



k"^ 



- Q 



Pa^^^i.) -T pp{kij_) 



A-2 



A-2 



(78) 



3.5 Summary of results for 

We can now collect our results for the components of the field A'^ and examine 
these more closely. By construction, this solution obeys d/^A^^ = 0. We have 
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also checked that it has the correct lowest order limit. It is useful to put this 
solution in the following compact form; 



+ 



+C^M [V{k,^) - {2^f5{k,^)] Y-^^ , (79) 



where we have defined the 4- vectors ^ as follows, 

C+{q,k^^) ^ , C-{q,k^^) ^ , Cj,(«,fei^) ^ -2kl , 

C+{q) = 2q+ , C-iq)^2^-2q- , = 2q\ (80) 

These coefficients are simply related to the well known Lipatov effective ver- 
tex [39-42] through the relation = + ^C^. For any momentum q, 
these 4-vectors satisfy 

q.C^=q.C^=0, (81) 

which ensures that the covariant gauge condition is trivially verified. Moreover, 
for a momentum q on-shell {2q^q~ = q'j_), we have the following additional 
properties: 

Cl = -4%M^ = d . (82) 

Remarkably, this ensures that when we use this solution in order to compute 
the production of on-shell gluons, we get only the correlator (UU'^) (the mixed 
correlators lUV'^) and (yC/^^, as well as IVV'^), do not appear). Moreover the 
non-zero scalar product is the square of the ordinary Lipatov vertex. This 
result suggests, as we will prove in the next section, that we likely have the 
property of fcj^-factorization for gluon production in pA. 



4 Classical gluon production in pA collisions 

4.1 fcx-factorization in gluon production 

From the results derived in the previous section, it is fairly easy now to derive 
the average multiplicity of gluons produced in pA collisions, to all orders in the 

source describing the nucleus. Using standard reduction formulas, the gluon 
production amplitude is given in the classical approximation by; 

Mx{.q) = q^A>^{q)e(^\q), (83) 
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where e^'^-' (q) is the polarization vector for a ghion of S-momentum q in the 
polarization state A. As usual, what this formula tells us is that we need to 
Fourier transform the gauge field, amputate the last free propagator (that's what 
the prefactor docs), put the result on shell, and project on the appropriate 
polarization vector. When we square this amplitude, we can sum over all four 
polarization states including the two non physical ones. Indeed, the property 
Qf^A^^ — ensures that the longitudinal polarization states are not going to 
contribute to the result. We can therefore use the identity, 

^e(^)(q)eW*(q) = -5^.. (84) 

X 

The average number of gluons produced in a collision is given by the relation 

^ J (27r)32i?, (2^)2 (27r)2 kl^k[\ 

X {Plaikl±)pp,a'{k[^)) {U\k21_)U{k'^^))^^, , (85) 

where (• • • ) denotes the average over the color sources. Note that the field A'^{q) 
of the proton alone does not contribute to on-shell gluon production. In order to 
rewrite this expression in a more intuitive way, let us first express the correlator 
of the proton sources Pp in terms of the unintegrated gluon distribution in the 
proton: 



5^(pla(fel-L)Pp.a'(fc'l±)) 



fcl_L + fe'ii 



(PXx_ 



■ (86) 



The integration over X ±^ runs over the transverse profile of the proton, = 
_ 1 is the dimension of the adjoint representation of SU{N), and dipp/ <P X ±_ 
is the number of gluons per unit area and unit of transverse momentum in the 
proton (i.e. the proton non-integrated gluon distribution per imit area). Note 
that the impact parameter does not appear here simply because we chose the 
origin of transverse coordinates at the center of the proton (6 will therefore ap- 
pear in the correlator involving the ?7's). This correlator can be approximated 
further if we realize that the support for the integration over X ±^ is the trans- 
verse area of the proton, a domain of typical size 1 / ^q^d ■ This means that 
the difference fei^ — k'^^, which is the conjugate variable of X ^, is at most 
of order A^^^ . Therefore, as long as we are interested in transverse momenta 
much larger than Kqcd, we can always approximate feij_ « k'l^. In particular, 
we can write: 
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Wc will also neglect the difference fcij^ — k'^^ in the other factors of eq. (85). 

Similarly, we can hide the correlator (U]^{k2±)U {k'2±)) in a function (^^ 
that describes the nucleus. By analogy with eq. (87), we define it as follows: 

<C/t(fc.,)f/(fc^,))„„, = ' 



X / d^Y^e^'^'^M-^. #.(fc.jy^-b) ^ ^gg^ 

where b is the impact parameter of the collision. The normalization of 9?^ has 
been chosen so that it is consistent with eq. (87) in the case of a dilute target. It 
is important to note that the function 97^ is not the canonical unintegrated gluon 
distribution of the nucleus, i.e. the expectation value of the number operator 
a^a^: it only coincides with the unintegrated gluon distribution at large fe_L. 
The correct interpretation of the function is that it is the square of the 
scattering amplitude of a gluon on the target nucleus. 

Neglecting the difference fci_L — A;']^j_ everywhere except in the exponentials, 
and using the explicit expression of , we obtain: 

^ 16a,7V ( d^q d^ki_ f,^^ dipp{k^\X ^) dip^iq^-k^\X ^-b) 



' T^d^ql J (27r)32£g (27r)2 J " d?Xi_ d^X± 

(89) 

This is the fc^-factorized form of the gluon multiplicity for pA collisions. We 
see that one can conveniently hide all the multiple scattering effects in a gener- 
alization of the definition of the unintegrated gluon distribution of the nucleus, 
the hard (pcrturbative) part of the matrix element remaining the same as in pp 
collisions. 

This expression for proton-nucleus collisions was first derived by Kovchcgov 
and Mueller [35], though not expressed in this form (see also Kopeliovich et 
al. [51] and Braun [46]). It was re-written in this form first by Kovchegov 
and Tuchin for deeply inelastic scattering [44] and subsequently by Kharzeev, 
Kovchegov and Tuchin for the pA case [45]. 



4.2 Compcirison with Dumitru and McLerran 

In this section, we compare our result for gluon production with the result 
obtained by Dumitru and McLerran [43] in the Fock-Schwinger gauge {x~^A~ + 
x~ = 0) which is an interpolation between two light-cone gauges. It is 
simpler to do this comparison at the level of the formula given by eq. (85). The 
latter may be rewritten as 

d^g ^ \_ fcPk^cfk[^ C„{q,ki^)-C^{q,k[^) 

d^q^dy Wtt^ J {2nf (2^)2 kf^K^ 

X (pL(fci±)Pp,a'(fe'i^)) {uHk2±)Uik',J)^^, . (90) 
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In the paper by Dumitru and McLerran, the gluon spectrum is given by 



dN„ 



(Pq^dy 



jtrc (|ai|^ + |a2p) , 



(91) 



with coefficients ai and a2 (converted to our notations) defined to be 

g3i7r/4 



p3i7r/4 



-5'iaf{XA_)di(f>U{XA_)ba) , 



a^iq^) = ^7^y d^x^e-'''^-^^e'^di(ar{x^)t'U{x^)ba) ■ (92) 

t'' is a generator of the fundamental representation of SU (N) and trc denotes a 
trace of color matrices. In these equations, ai^{x±) is defined by the relation 



(93) 



It is convenient to replace U{x±) and ai^{x±) by their Fourier transforms in 
the expressions for ai and a2- 



Uixi_) 



(27r)2 
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-C/(fe2±) , 



"1 (a;x)-y fc2^^f.«^'*= 



(94) 



We can thus write ai and 02 as 



g3i7i-/4 

„3i7r/4 



d^ki 

J2 



^<S'J-5Mpp,,(fcix)i''[/(fc2±)6 



/r2 



3i7r/4 [■ rfiU, , jUi p 

a2(gj = ^^™Pp,a(fei^)i''C/(fe2x)6a , (95) 

where now we have implicitly fe2± = Q^_ — ^i-L- Squaring these coefficients, 
taking the trace of their sum, averaging over the sources and inserting the result 
in eq. (91), we obtain: 



dNa 



d'^q±dy 



1 



hi 1.3 Uk Ul 
12 1 2 



U2 U2 



(27r)2 (27r)2 

X {pia{kii^)Pv..a'{KJ) {UHk,^)U(k',^))^^, . (96) 
At this point, it is trivial (although a bit tedious) to verify the identity 



-c^ [q, fei±) • {q, fe'i_L) = -2- 



S^jgkl ^ij^kl 



'>'l'>'2"'l "'2 ' 



(97) 



if the vector q is on-shell and fci + ^2 = fc^ + = g. This ends the proof that 

our result is equivalent to the result obtained by Dumitru and McLerran. Note 
that this equivalence works regardless of the model one chooses for the functional 
average over the hard sources. One need not specify this average to prove our 
result, nor even assume translational invariance in the transverse plane. 
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5 Cronin effect in classical gluon production 



The Cronin effect was discovered in proton-nucleus collisions in the late 70's [62- 
64] . The effect observed was a hardening of the transverse momentum spectrum 
in proton- nucleus collisions, relative to proton-proton collisions, that sets in at 
transverse momenta of order k± ^ 1 — 2 GeV, and disappears at much larger 
kj_'s. A corresponding depletion was seen at low transverse momenta, accom- 
panied by a softening of the spectrum. At that time, and indeed subsequently, 
the effect was interpreted as arising from the multiple scatterings of partons 
from the proton off partons from the nucleus [65]. As a result of such scatter- 
ings, the partons acquire a transverse momentum kick, shifting their momenta 
from lower to higher values, hereby causing the observed respective depletion 
and enhancement. At high k±, the higher twist effects, which, in the language 
of perturbative QCD, are responsible for multiple scattering [66,67] are sup- 
pressed by powers of k±. The relative enhancement of the cross-sections at 
moderate A:_l's should thus die away - and indeed, the data seemed to suggest 
as much. Though a qualitative understanding of the previously observed Cronin 
effect was suggested by perturbative QCD, a quantitative agreement for all its 
features (such as, for instance, the flavor dependence) is still lacking. 

The high energy deuteron-nucleus collisions at RHIC (and at the LHC in 
the near future) add another dimension to the Cronin effect, namely, its energy 
dependence. At high collision energies, the projectile probes the small x partons 
in the nuclear wave-function, whose distribution may be modified by the nuclear 
medium. As the energy increases, there are effectively more of them to scatter 
from and this should lead to an enhancement of the effect at higher energies. 
Eventually, the number of gluons saturates, implying a corresponding saturation 
of the effect of multiple scatterings. However, what happens also at high energies 
is a change in the spectrum and, as we shall see, that modifies qualitatively the 
net effect of the multiple scatterings. These various types of behaviors will be 
illustrated in this section. 

First data from RHIC provide indications on how the Cronin effect is mod- 
ified with energy or, cquivalently, with the rapidity. The x values probed in 
these experiments, at /c^ ~ 2 GeV, range from 10~^ in the central rapidity re- 
gion down to 10""^ at very forward rapidities. The most dramatic result is that 
of the BRAHMS experiment [1] which has taken data up to pseudo-rapidities 
7] = 3.2 . It is observed that the Cronin peak shrinks rapidly with rapidity 
and at higher rapidity, one sees that there is a significant suppression instead. 
Equally interesting is the centrality dependence of the effect [1]. At central 
rapidities, one observes that the Cronin peak is enhanced in more central colli- 
sions, while, for forward rapidities, the trend is reversed: more central collisions 
at forward rapidities show a greater suppression than peripheral collisions! 

The Cronin effect has been studied in the Color Glass Condensate frame- 
is The trends seen by BRAHMS are also corroborated by PHOBOS, PHENIX and STAR 
in more limited kinematic ranges [2-4], One can expect in the near future results from 
more detailed studies of the forward region in Dcuteron-Gold collisions at RHIC from all 
experiments-in particular from STAR and PHENIX. 
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work by a number of authors. These studies range from semi-quantitative [47, 
48,50,52-54] to detailed numerical solutions of the Balitsky-Kovchegov equa- 
tion [55]. There is a consensus among all these authors that the Cronin effect 
is suppressed at forward rapidities. It was shown early on by Dumitru. Jalilian- 
Marian, and one of us [47,48,50] that the McLerran-Venugopalan model shows 
a greater Cronin enhancement for larger values of the saturation scale Qs- In 
particular, the location and the width of the Cronin peak are roughly propor- 
tional to Qs. The saturation scale can be larger for greater centralities or for 
smaller values of x. The model therefore correctly predicts the greater enhance- 
ment for more central collisions for 77 = at RHIC. However, it fails badly 
when one goes to forward rapidities. Because x is small, the saturation scale Qg 
is larger and one gets more enhancement - rather than the observed suppres- 
sion. One possible cause for the failure of the McLerran-Venugopalan model 
is that it does not include the quantum evolution of the parton distributions. 
In fact, it was noticed early on that quantum corrections to the McLerran- 
Venugopalan model arc large [68,69]. The proper treatment of these corrections 
was performed by several authors - known collectively by their acronyms as 
JIMWLK [15-27], leading to the evolution equation for the weight functional 
W^[a;Q,/9^] for the sources p^. This equation can equivalently be re-cxprcssed 
as an infinite hierarchy of ordinary differential equations for n-point correlators: 
{UU) {UUU) , {UUU-- ■ ), where ?7 is a Wilson linens in the background field 
created by the source p^. For the 2-point fimction (^UW), a closed equation 
exists in the large limit: it is the Balitsky-Kovchegov equation [21,23]. This 
equation reduces to the well-known BFKL equation in the low density or large 
transverse momentum limit. Both the BK equation and the JIMWLK equation 
have been studied in analytical approximations (see [70] and [26] respectively). 
The BK equation has been studied numerically [71,72] and has been applied to a 
recent study of the Cronin effect [55]. The JIMWLK equation has been studied 
analytically in the "super saturated" small x region [73]. Detailed numerical 
solutions have only become available recently [74]. 

In this section, we will discuss the different features of the Cronin effect in 
the Color Glass Condensate framework. In all cases studied, the weight func- 
tional is a Gaussian. Then, the Cronin effect can be equivalently understood 
in terms of independent multiple scatterings, and the explicit connection be- 
tween the present formalism and Glauber's formalism of multiple scatterings is 
recalled in appendix C. We will first consider in detail the classical (no quantum 
evolution) McLerran-Venugopalan model's predictions for the Cronin effect. We 
next include quantum evolution "naively'' by introducing "by hand" an x de- 
pendence in the saturation scale: Qs Qaix), along the lines first suggested 
by Golec-Biernat and Wiisthoff based on their parameterizations of the HERA 
small-x data [75,76]. We show that this gives the wrong behavior for the Cronin 
effect at forward rapidities. We next consider analytical mean-field solutions of 
the JIMWLK equation - in a saturated regime where the system has no mem- 
ory of its initial conditions at small x. It was shown by lancu, Itakura and 

i^For more details, see [28-30]. 



26 



McLcrran [73] , that the weight functional of the sources in this regime is also a 
Gaussian, albeit a non-local Gaussian distribution. This Gaussian solution can 
be incorporated in our study of the Cronin effect, and we find that if one is in 
this "super saturated" regime already at y = 0, the Cronin effect is suppressed. 
Further evolution to y = 3 increases the suppression, but not by a great deal. 
We will discuss the implications of this study at the end of this section. 



5.1 Plain MV model 



If we consider a large nucleus, there is an approximate translation invariance of 
the correlator (UU'^^, which imphes that the function dipj^/d'^X±_ is approxi- 
mately independent of the location X_l in the transverse plane: 



A 7,2 



where we denote: 



(98) 



(99) 



This object is identical to the quantity C{k±) introduced in a paper by one of 
us and Peshier [49,77], except that it is now a correlator of Wilson lines in the 
adjoint representation. (This quantity, again for fundamental Wilson lines, was 
also studied in [14,78] where it was denoted 7(/c_l).) In the MV model, where 
the functional is a Gaussian weight given by 



J2^ PaAxi-)PaA^^) 
a x±_ — 



2/^ 



(100) 



this correlator can be computed in closed form. We get: 
(C/t(0)[/(a;i))^, = 5„6exp J d'y^Goiy^) - Go{y^ - x^) 

where Go is the 2-dimensional massless propagator: 

f d'^k± gi''±-(<^±-y±) 

Goi-±~y±)^J ^2 • 

In terms of the function C, the gluon spectrum reads: 



dN„ 



1 



d^q^dy IQ-n^q]^ J (27r) 



d^k 



-L l2 



kiC{ki_)ipp{q^ - fe_L) . 



(101) 

(102) 

(103) 



At large transverse momentum fc_L, we have the following behavior for C{k±): 



kiC(k^) 



(104) 
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which is the standard brcmsstrahhmg pcrturbativc taiL Note that this tail 
is proportional to /x^ oc A^^'^, where A is the atomic number of the nucleus. 
However, at lower transverse momentum (fc^ of the order oi ^ ~ g'^N^'^^ or 
smaller), the 1/fc^ growth stops and the fimction k'^C{k±) remains bounded 
when 0. This function was evaluated numerically in [49] and the typical 
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Figure 7: Behavior of (p{k±) oc k^C{k±) in the McLerran-Venugopalan model. 
In the upper right corner is plotted the quantity k\C{k±), which shows that 
the tail of the function (p{k±) behaves like In this plot, we have taken 

(5s_^ = = 25A^^^, where h-Qco is the infrared cutoff (which sets the 

momentum scale). 

behavior of k']_C{k±) is illustrated in figure 7. 

At low collision energy, we don't really know what the proton unintegrated 
gluon distribution is because it is mostly a non-pcrturbativc quantity. In order 
to perform a quick numerical study of eq. (103), we take the proton unintegrated 
gluon distribution ifip{k±) to be proportional^^ to k'j_C{k±), with the same func- 
tion C{k±) but now evaluated with a much smaller density /i^ (or saturation 
scale Ql^p)- We have now a C^{k±) and a Cp{k±), corresponding respectively 
to /U^ and /ip, and the ratio of the two saturation scales is taken to be: 

4 = « 6 . (105) 

^^Obviously, using a proton gluon distribution proportional to k^C{k±), i.e. to a correlator 
of Wilson lines, is a totally unjustified ansatz. Indeed, we have only proven fcj_ -factorization 
for a calculation performed at the lowest order in the source describing the proton. However, 
this ansatz for ipp has the expected large kj_ tail, and here one could simply see it as a way 
of introducing an infrared cutoff (necessary because the bremsstrahlung spectrum in l/fc^. 
not integrable at low k±). Moreover, what we take for (pp is not essential for the discussion of 
the Cronin effect, as it is mostly due to properties of the small-o; wave function of the target. 
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Up to a constant normalization factor, we have 



dN„ 



(X 

pA 11 



<Pq±dy 

for proton-nucleus collisions, and 

dN, 



(106) 



d?q^dy 



oc 

pp 11 



^ / ^^^^ ~ C,{k^)CM^ - k±) , (107) 



for pp collisions. The spectrum in pA collisions can be estimated numerically 
using eq. (106), and we have displayed the result in figure 8. One can see clearly 
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Figure 8: Numerical calculation of the spectrum in pA collisions. We use 
Qs,a/^qcd = 5 and Qs,p/^qcd = 2. The dotted curves are fits to power 
laws. 



on this plot a change of behavior between a tail in l/q\ and an intermediate 



region where the slope of the spectrum is l/q\, as predicted by Dumitru and 



McLerran. 

We can also estimate the "nuclear enhancement ratio" RpA as follows: 



dN„ 



Neon d,^q±d.y 

RpA = 



1 lSw{^±-k±mCAq^-k^)Cp{k^) 



pA _ J- J (27r) 



I ^ili. - k±ni Cp{q^ - k^)Cp{k^) 



(108) 

We have computed numerically the ratio of these integrals, and the result is 
plotted in figure 9. One can see clearly the Cronin effect as a small bump at 
q±/^Qcn ~ 11- ratio is smaller than unity at small momenta, and then 
remains larger than 1 at larger momenta. 
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Figure 9: The ratio RpA in the McLerran-Venugopalan model, without quantum 
evolution. We use the values Qs,a/ ^qcd = '^•9 ^^id Qs,p/^qcd = 2, which 
corresponds to a ratio /i^/yUp « 6. 



When the saturation scale of the miclcus is much larger than that of the 
proton, the location of the Cronin peak can be estimated in a semi-quantitative 
way as follows. A good approximation of the behavior of the function C^{k±) 
for transverse momenta k± comparable to the saturation scale or smaller is: 

C^^k^)^^e-'-'-/<^, (109) 

with = QlJniQlJAl^J and Ql^ = Sn^a^.Nnl. Assuming that the 
ratio RpA can be approximated by the ratio of the distributions in the targets^^ 

and that Cp{k±) cx l/kj_ in the relevant region^^, we see that RpA has a max- 
imum at k± = Qs^/a/tt- The apparent mismatch with the location of the 
maximum in figure 9 is due to an insufficient separation between Qs,a 9'nd Qs,p 
(in addition to the fact that the two saturation scales are not very large com- 
pared to A^^^). 

^®This would be exact if the unintegrated gluon distribution of the projectile were propor- 
tional to a i5(fex)- In practice, this is a reasonable approximation if the saturation scale in 
the nucleus target is much larger than that of the projectile. 

^^In the regime where the saturation scale in the proton is much smaller than that in the 
nucleus, the bulk of the gluon distribution of the nucleus corresponds to momenta that axe in 
the tail of the gluon distribution of the proton. 



30 



5.2 Quantum evolution introduced via //^(x) 

The results of the previous subsection have been obtained in the plain McLerran- 
Venugopalan model, which does not contain any quantum evolution. This means 
that the quantities /x^ , /x^ (i.e. the saturation scales of the nucleus and of the 
proton) do not depend on the momentum fraction x. One consequence is that 
the gluon multiplicity is exactly boost invariant - independent of y. A naive ex- 
tension of the previous calculations is to assume that the only effect of quantum 
evolution is to change the values of /i^ and /i^ by making them x dependent, 
while preserving the local Gaussian nature of the functional weight. Following a 
popular parameterization, originally due to Golec-Biernat and Wiisthoff [75,76], 
we parametrize the saturation scale as 

QlJQl = A"'{^)\ (111) 

with Qq = 1 GeV. We use the values of Xq and A obtained in the recent paper 
by lancu, Itakura and Mimier: xq = 0.67 • 10~^ and A = 0.25 from their fits to 
the HERA F2 data [79]. The values of x in the proton and in the nucleus are 
given by: 

r - ^P+V 

X, = ^^e-y , (112) 

where ^/s = 200 GeV is the energy per nucleon in the collision. For definiteness, 
we set the infrared cutoff to A^^.^^ = 200 MeV. We then simply repeat the 
previous calculation with the following formula for RpA- 

I$^{q±-k^rklCp{x,,k^)Cp{xp,q^-k^) ' 

The results of the numerical evaluation of this ratio are displayed in figure 10. 

The characteristic Cronin bump is still present both at y = and at y — 3. 
The primary difference is that at larger rapidity the maximum is more pro- 
nounced and is pushed towards larger momenta. This is a consequence of the 
fact that it is the Qg of the nucleus that determines the location of the maximum 
to a large extent (see eqs. (109) and (110)). Similarly, more central collisions 
would give a greater enhancement as well, since they correspond to a larger Qg. 
This picture is consistent with the centrality dependence observed at y = in 
Deuteron-Gold collisions at RHIC. It is however in complete disagreement, at 
forward rapidities, with the observation made recently by the BRAHMS exper- 
iment [1]. 
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Figure 10: The ratio RpA in the McLerran-Venugopalan model, with quantum 
evolution introduced via the x dependence of the saturation scale. 



5.3 Non-local Gaussian distribution 



One can generalize the previous model for the correlator (j7^(0)C/(a;^))^j^ by 
allowing transverse non-localities in the Gaussian functional Wj^ . This amounts 
to generalizing eq. (100) by allowing /u^ to depend not only on x but also on 
the transverse coordinates: 



^aI^^Pa] = exp 



J 2nj^{x,x±-yj_) 



so that the average of a product of two p is given by 

{PA,a{xi-)PA,biy±)) = SabP'lix,X± - yj_) . 



(114) 



(115) 



In this model, we have assumed that the nucleus is invariant by translation in 
the transverse direction, which implies that /^^ depends only on the difference 
x± — y ^. For instance, this extension can be used to incorporate effects of 
color charge neutralization [80] by making the function ii^^{x,x± — y_|_) vanish 
at distances larger than a certain scale. In this model, one can still calculate in 
closed form the correlator of two Wilson lines: 



{U\0)U{x^))^^ = Sabe^p 



I N 

— I d'^y_i_(fz_Ltil{x,y_i_-z±) 



: (Go(yi) - Go{y^ - x^)) {Go{z^) - Go{zi_ - x^))\ , 

(116) 
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or by going to Fourier space, 



2-K 



+ OC 

i'- 



dk± 



(1 - Jo{kjLXA_))ii\{x,ki_) 



with 



{x,k±) 



I 



(117) 



(118) 



(We have assumed rotational invariance in the transverse plane. This Fourier 
transform therefore depends only on fcj_ = |fc_L|). 

Such a non-local Gaussian distribution of color sources has in fact been 
predicted by lancu, Itakura, McLerran [26] , as a mean- field asymptotic solution 
of the JIMWLK evolution equation for the deeply saturated regime. In this 
asymptotic regime, the color neutralization is controlled by saturation physics 
and therefore happens at the scale Q 
given by: 

" ^-ii\{x,ks_) = —k'i In 



instead of A„p^ , and the function fj? is 



27r 



1 



kl 



(119) 



In this equation, 7 is some anomalous dimension (7 k, 0.64 for BFKL) and 
c « 4.84. Before going on, it is instructive to discuss the unintegrated gluon 
distribution ip{x,k±_) oc k^C{x,k±) for this model. The function C{x,k±) can 
be expanded in powers of /i^(x, fcj^), which corresponds to an expansion in the 
number of collisions (see the appendix C): 



C{x, kj_) = Ci{x, kj_) + C2{x, kj_) + ■ 



(120) 



ith^o 



C.ix, k^)=9 = Tcl^J^y+[-\ 




(121) 



and 



C2(a;,fcj_) 



d Pi_ 
(27r)2 



o^a(^'^-l) y'\{x,Pl.) 



ki 



Pi 



(122) 

One can note that with the /i^(.T,p^) given in eq. (119), this quantity is in- 
frared finite. A numerical study of this 2-scattering correction shows that it 
is negative at low k±_ and positive at large k_\_. Therefore, the property that 
rescattcrings move particles from low k±_ to high k± regions of the spectrum 
holds also in this model. However, since in this model the spectrum has a 
smaller slope at large k±_ {k']^^^'^'^^ with 7 w 0.64) than the MV model spec- 
trum (fcj^), this effect of rescatterings is much less important. That higher 

^"Note that tp{x, k±) depends on x and k±_ only via the ratio k±/Qs{x). 
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Figure 11: The function k'j_C{x, k±) in the model defined by eq. (119) (plotted 
against k±/Qs,A{x) thanks to its scaling properties). The insert in the upper 
right corner represents the ratio C(a;, A;_L)/Ci(a;, fcj,). 



twists tend to push the spectrum towards higher momenta is clearly seen in 
figure 11, where we have both represented the function k'^C{x,k±) (in the main 
plot), and the ratio C(x, k±)/Ci{x, k±) (in the smaller insert). However, as the 
plot of C{x, k±)/Ci{x, k±) also shows, the maximum over the single scattering 
term is now of only 3%, compared to 75% in the MV model (see the figure 7). 

Assuming that the nuclear enhancement ratio RpA is well approximated by 
eq. (110), we can obtain its asymptotic value at high k± from the leading term 
Ci(a;, fc_L): 

RpA « A-'^'(%A =A(''-'^/K (123) 

Since 7 < 1, this asymptotic value is smaller than 1. For 7 « 0.64 and ^ ~ 6, 
we expect RpA ~ 0.52 at large momentum (this ratio would probably go to 
1 asymptotically if we let the anomalous dimension 7 depend on transverse 
momcntiim and go to 1 at large momentum). Moreover, the ratio RpA in this 
approximation is the ratio of two functions similar to the one plotted in the 
upper right corner of figure 11, rescaled horizontally by their respective Qs and 
vertically by Q^"* ■ Therefore, we expect the ratio RpA to be below its asymptotic 
value A^'''"^)/^ at small k±, to eventually become larger than A^''~^^^'^ and reach 
a very mild maximum (of a few percent), and finally to reach the asymptotic 
value from above. 

We have studied the ratio RpA numerically, and the rcsidts are displayed 
in the figure 12. The asymptotic value of RpA is in relatively good agreement 
with the crude estimate given in eq. (123). It seems also that the very mild 
maximum we predicted on the basis of the approximate formula of eq. (110) has 
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Figure 12: The ratio RpA in the non-local Gaussian model of eq. (119). 

been washed out by the convolution present in the exact formula. The shape of 
the curve now agrees qualitatively with BRAHMS data for y = 3. However, it 
also gives too much of a suppression at t/ = 0! 

5.4 Interpretation of results 

We can summarize the key results of the model calculations in this section as 
follows. 

• Tree level partonic re-scattcring responsible for the classical Cronin effect 
in high energy deuteron-Gold collisions is computable in the color glass 
condensate approach and when the weight for the color sources is Gaus- 
sian, the formalism reduces to the familiar Glauber formalism of indepen- 
dent multiple scatterings. In the simple McLerran-Venugopalan model, a 
peak appears at k± ~ Qs,a^ which is more pronounced for more central 
collisions. 

• Quantum corrections (due to small-a; evolution) to the tree level partonic 

re-scattering can be "naively" included by letting Qs Qs{x). The pa- 
rameterization of Qs {x) as a function of x that we have used is of the form 
first suggested by Golec-Biernat and Wiisthoff but with the parameters 
from a recent fit by lancu, Itakura and Munier to the HERA data. We 
call this model of quantum evolution naive because it does not consider 
the possibility that quantum evolution will change the form of the weight 
functional from the Gaussian form of the MV model. This naive quan- 
tum evolution model predicts a larger Cronin effect at larger rapidities 
(because Qs grows as x decreases) and is in sharp disagreement with the 
RHIC d-Au data at forward rapidities [1]. 



35 



• Proper quantum evolution in the CGC is described by the JIMWLK evo- 
lution equation. For the two-point correlator, of interest in Cronin studies, 
the large N limit is described by the Balitsky-Kovchegov equation. Like 
any other evolution equation, these equations evolve from an initial condi- 
tion at some large x = xq. At very small a; <C xq, the system loses memory 
of its initial conditions. The JIMWLK equation, in this limit, is well de- 
scribed by a mean field approximation, and analytic expressions have been 
derived. Interestingly, the weight functional in this extreme limit is also 
a Gaussian, but with a non-local kernel. That is, the variance /u^ of the 
Gaussian sources is now a function of k±: /i^ (x) fi'^^ {x, k±). At low k±, 
A*A ^1. ^^(^ + (Q^/^J-V)^ where 7 is an anomalous dimension (7 0.64 
for BFKL evolution). Because the weight is Gaussian, the physics remains 
that of multiple independent scatterings. However, as our calculation in 
section 5.3 shows, the effects of multiple scatterings depend very much 
on the shape of the spectrum. It follows that the qualitative behavior of 
the ratio RpA is different from what it is in the MV model with a local 
weight. If we assume that y = and y = 3 at RHIC lie in this "super 
saturated" regime, we find a suppression both at y = and y = 3. Thus 
while our results for y = 3 are qualitatively close to the BRAHMS data, 
the y = results are decidedly not close to the corresponding y = data. 
Our result therefore suggests that quantum evolution is not yet significant 
at J/ = 0. This of course was also evident from the disagreement of the 
prediction of Kharzeev, Levin and McLerran with the RHIC d-Au data 
at y = 0. Our result merely provides additional insight into why one has 
this disagreement. 

From the comparison of the results of our simple studies with the qualitative 
features of the RHIC d-Au data emerges the following picture, which will have 
to be substantiated by more detailed (read quantitative) computations and a 
thorough confrontation to the experimental data. That the MV model works 
reasonably well at x ~ 10~^ - central rapidities at RHIC - and reproduces 
key qualitative features of central d-Au collisions suggests that it has the right 
physics built in and that it could be a good model of the initial conditions 
for quantum evolution as a function of rapidity. The fact that at the same 
time the MV model fails badly at forward rapidities is suggestive of important 
(and non-trivial) quantum evolution effects that cannot be accounted for simply 
by a rescaling of the saturation momentum. This interpretation is consistent 
with the fact that the asymptotic mean-field solution of the JIMWLK evolution 
equation is in qualitative agreement with the d-Au data at j/ = 3. In particular, 
a prediction of this mean-field solution is that the ratio RpA tends to A^'^~^^/^ 
at large transverse momentum, where 7 is the BFKL anomalous dimension. 
If there is some truth to this interpretation, one would also predict that the 
suppression of the ratio RpA would change very little if even larger rapidities 
were available (since y = 3 seems to be already in the asymptotic regime) . This 
corroborates the numerical results obtained previously from solutions of the BK 
equation by Albacete et al [55], and is also confirmed by lancu, Itakura and 
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Triantafyllopoulos [81]. The preliminary RHIC d-Au data appear to similarly 
show a rapid change at lower rapidities but a weaker one at the higher rapidities 
probed. 

6 Summary and Outlook 

We have performed in this paper a systematic study of gluon production in pA 
collisions within the framework of the Color Glass Condensate. Working in the 
covariant gauge, we have solved the Yang-Mills equations to lowest order in the 
dilute proton source but to all orders in the nuclear source. Our explicit results 
for the gauge fields arc useful in computing quark production in pA collisions. 
Our results on these are reported on in a companion paper. We have computed 
gluon production and demonstrated explicitly how fcj_-factorization arises in this 
approach. We have demonstrated explicitly the equivalence of calculations of 
gluon production in different gauges. 

The fc_L-dependent distribution that describes the nucleus is proportional to 
a two-point correlator of Wilson lines. This suggests how small-x quantum evo- 
lution effects can be taken into account in studying proton-nucleus collisions in 
the CGC framework since the two point correlators satisfy non-linear evolution 
equations with x, or cquivalently the rapidity y. We next studied in some detail 
the Cronin effect in the CGC framework. A rapid transition from the regime 
where multiple scattering effects are dominant to one where quantum evolution 
effects arc dominant is suggested by simple analytical studies corroborating pre- 
vious numerical studies. If confirmed by additional experimental results, this 
rapid transition has important ramifications in the near future for heavy ion 
experiments at the LHC and even for the proposed electron-heavy ion deeply 
inelastic scattering experiments (eRHIC) proposed at BNL. 
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A Proof of equation (26) 

Let's assume we want to solve the equation: 

+ iF{x)d+)A{x) = J{x) , (124) 

with a source J{x). The function F{x) is some Hermitian matrix- valued object, 
and we assume that it does not depend on x~, i.e. that it commutes with the 
derivative 9+. We want to find the function A{x) for times a;+ > 0, given 
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the value of A{x^ = 0,x~,x±) (and maybe of its first derivatives, since the 
equation is of second order in derivatives) and the value of A(x) and its first 
derivatives on the boundary of the "spatial" directions (here x~ and x±). In 
order to do that, we must start with the equations for A{x) and for the free 
retarded propagator^ ^: 

iOy +iF{y) d+)A{y) = J{y) , 

Gn{x,y){ny -iF{y) '^+) = 5{x-y) . (125) 

We have written the equation for the propagator with respect to its second 
point, which imposes to put the term F{x) on the right because of its matrix 
nature. The arrows on the differential operators indicate on which side the 
derivatives act. Then, multiply the first equation by G^{x,y) (on the left), the 
second equation by A{y) (on the right), subtract them, and integrate over y 
(with y+ starting at 0). This gives: 

A{x) = J d^yG^{x,y)J{y) 

y+>0 

+ 

y+>0 

-G^{x,y){ny +iF{y) d+)A{y)} . (126) 



I d^y{G^ix,y)iny-iF{y)d+)A{y) 



This is a generalization to the case of eq. (124) of the standard Green's theorem 
of electrostatics. The first term is the one we would naively expect if the initial 
and boundary conditions for the field and its derivative are zero. In order to 
deal with the second term, we first write Uy = "^dydy — V^^^. For the term 
involving the transverse derivatives, we can write: 



^2 ^2 

GH{3^,y) A{y)-G^{x,y) A{y) 



v± 

G^{x,y) Vy^ A{y) - G^{x,y) Vy^ A{y) 



• (127) 



We then perform the integration over yj^ using Stokes theorem, and the corre- 
sponding term becomes: 



j dy^dy j> 
y+>o am^ 



dn_i_ 



Gn{x,y) Vy^ A{y) - G^{x,y) Vy^ A{y) 



(128) 



^^The equation for the propagator when one acts on the right has been written by analogy 
with what happens when the "external field" is constant: in this case, the propagator is 
invajriant by translation and depends only on x — y. Hence, the first derivatives with respect 
to y enter with the opposite sign compared to the first derivatives with respect to x, while the 
second derivatives are the same. 
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where 9]R^ is the boundary in the direction of transverse coordinates, and n± 
a unit normal vector on this boundary. Similarly, we can deal with the term in 
dyd~ by first writing: 

GA^,y) d+ dyAiy) - G, {x,y) ~d+ ~d- A{y) 

= d+{G,(x,y) 1)-A{y))- d-{G^{x,y) 3+^(2/)) • (129) 
Therefore, the corresponding term gives: 

3/~=+oo 
y- = -oo 



J dy+d' 



G^{x,y) dyA{y) 



y+>o 



-2 Jdy-d^y^_ [G^{x,y) d+A{y) 



y+=0 



(130) 



Using the fact that G^ is a retarded propagator, the boundary term at y+ = +00 
is zero. Finally, we need to deal with the term in F{y)d^ . We have: 

j d^y{G,{x,y)F{y) d+ A{y) + G ^{x,y)F{y) ^+ A{y)] 

I j dy+d^yAG,{x,y)F{y)A{y)]i:ztZ ■ (1^1) 



!/+>0 



y+>0 



Here, we have used the fact that F{y) commutes with 9+ and the fact that 
G^{x, y) is a retarded propagator. If we combine all the pieces together, we get: 



A{x) = j d^yG^{x,y)J{y) 

dy-d^y^ G^{x,y)2dtA{y) 



y+>o 



y+=o 



+ 



+ 



j dy+dy j> 

j/+>0 dm? 



dn± 



G^{x,y) Vy^ A{y) - G^{x,y) Vy^ A{y) 



dy+d'y^ G^{x,y){2 dy - iF{y))A{y) 



y =+00 



y =-00 



(132) 



y+>0 



If we assume that the field and its derivatives vanish fast enough at infinity in 
the spatial directions, we get the simpler formula: 

A{x) = I d^yG^{x,y)J{y) 

y+>0 

J dy-d^y^G^{x,y)2^+A{y) . (133) 

2/+=0 



+ 
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Eq. (26) is a special case of this formula, where the source term J is zero. 
Naturally, in the above derivation, the hyperplane = could have been 
replaced by any plane of constant y+ on which we decide to set the initial 
condition. 



B Simplification of equation (71) 

In order to simpliiy eq. (71), we need to compute the following convolution 
product: 



U (E>{A- ■T)(g)V = / dx+U{e,x+){A-{x+) ■ T)Vix+ ,0) . 



(134) 



We have not written the transverse coordinates in this expression; they are the 
same for the three factors and are irrelevant for the following discussion. Let us 
first define 



W{b+,a+\x+) =V+ 



expi- 





J dz+{l 



x+j)A-{z+)-T 



This object interpolates between [/(6+,a+) and V^(6+,a+): 

W{b+, a+|x+) = V{b+,a+) if x+ > 6+ , 
W{b+,a+\x+) = U{b+,a+) if x+ < a+ . 



We then have 

dW{b+,a+\x+) 
dx+ 



(135) 



(136) 



= -i^Wib+,x+\x+)iA-{x+)-T)Wix+,a+\x+) 
= -i^U{b+,x+){A-{x+)-T)V{x+,a+). 



(137) 



Then, integrating over a;+ from to e, we get: 



i^[U ® {A- ■T)®V]= U{e, 0) - V{e, 0) 



(138) 



Let us now go back to cq. (71). Using the fact that < i;"*" < e, we can 
approximate exp(ig"x+) w 1. Using cq. (138), we have at this point: 

q''A\q) = 2 jdv-d^v^e'^''^''~-'i^-''^\V{vA_)-U{vA_))diA+{v) 

+ ig jdv-d'^v±e''^i^''~-''^-''^^{[U (g)idlA- ■T)(E)V]iv^) 



-[U^{A--T)^{dlV)]{v^)}A;{v). 



(139) 
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The next step is to make more explicit the derivative dlV: 

diV = i^[V^{diA-.T)^V], (140) 

which leads to^^: 

[U {diA- • T) ® y] - [f/ ® {A- . T) ® {diV)] 

= [u^ {diA- • T) F] - i|[c/ ® (A; • T) ® y o {diA- ■t)^v] 

= ii[V<»{dlA--T)<S)V] = yiV. (143) 
z ig 



C Glauber interpretation of C{k±) 

It is fairly easy to write the Fourier transform 5acC{pi_) of (i7^(0)J7(a;j_)) in 
a form that makes its interpretation more intuitive in terms of the Glauber 
picture. Our starting point is cq. (117) (which is the most general form one can 
obtain in the case of a Gaussian distribution of ) : 

C{p^) = J d^x^e'P^-^^e^p J !^{l-Jo{k^x^))i,l{x,k^) 

■(144) 

At this point, it is convenient to factor out of /i^ a constant /Xg that sets the 
scale: 

^il{x,k±)=^ilf{x,kjL). (145) 

Typically, /Iq is given by = pL, where p is the density in the target and 
L is the longitudinal size of the target. Then, we can introduce the following 
shorthands: 



+ 00 



ftot 



/ ^f{x,ki_), 



27r J 



+ OC 



a{x^)=^ J ^Mk±x^)f{x,k^) . (146) 

^ 



We use the fact that the convolution product defined above is associative in the following 
sense: 

[F [G B H]] = [[F (8) A(^G] ^ B ^ H] . (141) 

This can be proven simply by permuting the integrals over the two intermediate coordinates: 

/ dx+ f dy+ = f dy+ f dx+ . (142) 

Jo Jo Jo Jy+ 
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(Ttot is the total cross-section of an incoming parton on a parton of the target, 
and a{x±) is the Fourier transform of the differential cross-section (with respect 
to the transverse momentum transfer). We can therefore write: 

/+00 2n 
(fx±e'P^-''^^\[a{x±)r . (147) 



If we now replace each factor of (j{x±) by its explicit expression as a Fourier 
transform of a{k±) = g'^N f(x, k±)/kj_, we get: 



+°° ,,2n r rl'^U. . ^2f^ 



2 XXr ,,2" f rl^lc 



n=0 
2s 



r)2 (27r)2 

X (27r)25(fei_L + • • • + fe„± - p_L)a(fei_L) • • • a{kn±) ■ (148) 
Finally, if we notice that: 

t^=p^ J dzij dZ2--- J dZn, (149) 

Zl Zn — 1 

we see that the object C{p±) can be interpreted as the Glauber form of the 
cross-section of a parton undergoing multiple independent scatterings in the 
target (the index n in the sum is the number of such scatterings, and the Zi 
are the longitudinal coordinates of the scattering centers) . Note that this corre- 
spondence is certainly only valid when the functional distribution of the source 
is a Gaussian distribution (having a distribution with non-zero (p^PaPa)' 
{PaPaPaPa} ' ' ' 1 introduces correlations between the different collisions and the 
Glauber picture of independent successive collisions breaks down). The form 
given in eq. (148) for C{p±) also makes obvious the following sum-rule: 

^Cip^) = 1 • (150) 

Indeed, integrating eq. (148) simply removes the delta function that was con- 
straining the variables ki±, so that integrating over each ki± now gives crtot- 
One can then rc-exponcntiatc the sum over n, which gives a factor exp(/ioO'tot) 
that compensates the prefactor exp(— /Uodtot)- 
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